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A reaction-diffusion model for the growth of avascular tumor 



O 
O 

> 
O 

o 

(N 



S. C. Ferreira Jumori,Q M. L. Martins^,! and M. J. Vilela^ 

^ Departamento de Fisica, Instituto de Ciencias Exatas, 
Umversidade Federal de Minas Gerais, CP 702, 30161-970, Belo Horizonte, MG, Brazil 
^ Departamento de Fisica, Universidade Federal de Vigosa, 36571-000, Vigosa, MG, Brazil and 
^ Departamento de Biologia Animal, Universidade Federal de Vigosa, 36571-000, Vigosa, MG, Brazil 

(Dated: February 1, 2008) 

A nutrient-limited model for avascular cancer growth including cell proliferation, motility and 
death is presented. The model qualitatively reproduces commonly observed morphologies for pri- 
mary tumors, and the simulated patterns are characterized by its gyration radius, total number 
of cancer cells, and number of cells on tumor periphery. These very distinct morphological pat- 
terns follow Gompertz growth curves, but exhibit different scaling laws for their surfaces. Also, the 
simulated tumors incorporate a spatial structure composed of a central necrotic core, an inner rim 
of quiescent cells and a narrow outer shell of proliferating cells in agreement with biological data. 
Finally, our results indicate that the competition for nutrients among normal and cancer cells may 
be a determinant factor in generating papillary tumor morphology. 

PACS numbers: 87.10.-(-e, 87.18.Hf, 87.15.Vv 



I 

o 
o 



> 

ON 

O 



X3 
O 

o 



X 



I. INTRODUCTION 



Cancer is a disease derived, with few exceptions, from 
mutations on single somatic cells that disregard the nor- 
mal controls of proliferation, invade adjacent normal tis- 
sues and give rise to secondary tumors (metastasis) on 
sites different from its primary origin §. Although can- 
cers are extremely diverse and heterogeneous, a small 
number of pivotal steps associated to both deregulated 
cell proliferation and suppressed cell death is required 
for the development of any and all tumors. Indeed, all 
neoplasms evolve accordingly to an universal scheme of 
progression |j . Neoplastic cells accumulate a series of 
genetic or epigenetic changes along the tumor progres- 
sion in response to natural selection and as an integrated 
defense program against stress situations similar to the 
response of bacterial colony facing severe and sustained 
threats jj]. But, unlike of the bacterial growth, tumor 
progression involves a complex network of interactions 
among cancer cells and its host microenvironment 
It is well known that stressed bacterial colonies can de- 
velop diffusion-limited fractal patterns |^. Recently, 
normal and tumor cell patterns in vivo and in vitro were 
characterised by their fractal dimensions and cluster size 
distribution functions [|[ ^, ^ , reinforcing the great cur- 
rent interest in the search for basic principles of growth in 
living organisms, which are the most complex and chal- 
lenging self-organized systems. In particular, for cancer 
growth, one of the most aggressive phenomena in biol- 
ogy, numerous mathematical models have been recently 
investigated. Examples include studies, based on classi- 
cal reaction-diffusion equations, of the growth of tumor 
spheroids [Gdl, cancer evolution and its interation with 




and the fundamental problem 
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the immune system 11 
of tumor angiogenesis 

Recently, we |l^ proposed a diffusion-limited 

model to simulate the growth of carcinoma "in situ" in 
which cell proliferation, motility and death are locally 
regulated by the concentration of growth factors pro- 
duced by each cancer cell. This model was able to gener- 
ate compact, connected and disconnected morphologies 
which progress in time according to Gompertz growth 
curves, and for which the tumor gyration radii scale as 
in the Eden model for an asymptotically large number 
of cells. These features were shown to be independent 
of the detailed functional form of the microscopic growth 
rules. In contrast, the structure of the tumor border is 
influenced by the nature of the growth rules as indicated 
by the different scaling laws for the number of periph- 
eral cancer cells. In particular, for disconnected and con- 
nected patterns the surface widths scale with exponents 
similar to those observed in bacterial colonies. Although 
the similarities between the simulated and histological 
tumor patterns were encouraging, the model was unable 
to generate papillary and ramified morphologies found in 
many of epithelial cancers and trichoblastomas. In con- 
sequence, we were lead to investigate the role of nutrient 
competition in neoplasic development, a biological fea- 
ture excluded in our previous model but central in the 
growth of fractal bacterial colonies. Indeed, cancer cells 
subvert the evolutionary adaptations to multicellularity 
and revert to a largely nutrient-limited style of growth. 
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In this paper we analyze the avascular cancer growth in 
a model including cell proliferation, motility and death as 
well as competition for nutrients among normal and can- 
cer cells. The cell actions (division, migration and death) 
are locally controlled by the nutrient concentration field. 
In Section 2 we introduce the cancer growth model. In 
Section 3 the simulational results for the tumor patterns, 
growth curves and scaling laws are discussed. In Section 
4 the diffusion of growth factors secreted by the cancer 
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cells are added to the nutrient-limited model. Finally, we 
draw some conclusions in Section 5. 



II. THE CANCER GROWTH MODEL 

The basic biological principles included in the model 
are cell proliferation, motility and death and competition 
for nutrients among normal and cancer cells. Nutrients 
(oxigcn, amino acids, glucose, metal ions, etc.) diffuse 
from a capillary vessel through the tissue towards the 
individuals (normal and cancer) cells. Under restricted 
nutrient supply the growth rate of cancer cells is limited 
by its ability to compete for nutrients with the normal 
cells. In our model the division, migration and death of 
each cancer cell is controlled by the nutrient concentra- 
tion in its local microenvironment. 



A. The tissue 

The studied system consists of a tissue fed by a single 
capillary vessel. The tissue is represented by a square 
lattice of size (L -|- 1) x (L -I- 1) and lattice constant A. 
The capillary vessel, localized at the top of the lattice 
at X = 0, is the unique source from which nutrients dif- 
fuse through the tissue towards the individuals cells. Al- 
though a tumor mass is composed of different cell sub- 
populations we shall consider only three types: nor- 
mal, cancer and tumor necrotic cells. Any site, with 
coordinates x = (iA,jA), i,j ~ 0, 1,2,... ,L, is oc- 
cupied by only one of these cell types. In contrast to 
normal cells, one or more cancer cells can pile up in a 
given site. In turn, tumor necrotic cells are inert and, 
for simplicity, will be considered always as a single dead 
cell. Thus, each lattice site can be thought of as a group 
of actual cells in which the normal, necrotic and can- 
cer cell populations assume one of the possible values 
(T„(x, i) = ad{x,t) = 0,1 and ac(x,t) = 0,1,2,..., re- 
spectively. As initial "seed" a single cancer cell in the 
half of the lattice {x = LA/2) and at a distance Y from 
the capillary vessel is introduced in the normal tissue, 
in agreement with the theory of the clonal origin of can- 
cer |l7| . Periodic boundary conditions along the horizon- 
tal axis are used. The row i = represents a capillary 
vessel and the sites with i = L + 1 constitute the external 
border of the tissue. 



The nutrients 



non-essential for cell proliferation, described by the con- 
centration fields N{x,t) and M{x,t), respectively. How- 
ever, it is assumed that both nutrient types have the same 
diffusion coefficients and consumption rates by the nor- 
mal cells. These concentration fields obey the diffusion 
equations: 



dN{x, t) 
dt 



DV'^N{x, t) - jN{x, t)an{x, t) 

-\NlN{x,t)<T,{x,t) (1) 



and 



dM{x, t) 
di 



^ DV'^M{x, t) -jMix, t)crn{x,t) 

-\MlM{x,t)n,{x,t) (2) 



in which the nutrient absorption terms are proportional 
to the cell populations present in each site, and it is as- 
sumed differentiated nutrient consumption rates for nor- 
mal and cancer cells by factors \n and Am- It is impor- 
tant to notice that the present model assumes the sim- 
plest form for the nutrient diffusion phenomena, i. e., lin- 
ear equations with constant coefficients. Also, \n > Am 
is used, reflecting the larger cancer cells affinity for es- 
sential nutrients. 

The boundary conditions satisfied by the nutrient con- 
centration fields are N{x = 0) = M{x = 0) = K^^ rep- 
resenting the continuous and fixed supply of nutrients 
provided by the capillary vessel; N{y = 0) = N{y = 
LA) and M{y = 0) = M{y = LA), corresponding to 
the periodic boundary conditions along the x-axis; at 
last, Neumann boundary conditions, dN{x — LA)/dy = 
dM[x ~ LA)/dy = 0, are imposed to the border of the 
tissue. The hypothesis that a blood vessel provides a 
fixed nutrient supply to the cells in a tissue is a simplifica- 
tion, which neglects the complex response of the vascular 
system to metabolic changes of cell behavior p9| . 

In order to reduce the number of parameters in equa- 
tions (1) and (2), the new dimensionless variables are 
defined: 



A2 ' A ' 
M' = ^,a = A,, 



(3) 



Using these new variables in Eqs. (1) and (2) and omit- 
ting the primes we obtain 



As considered by Scalerandi et al. |l8| , we assume that 
dividing cancer cells are especially vulnerable to some 
critical nutrients such as iron, essential for DNA syn- 
thesis and, therefore, for cell division. The many other 
nutrients necessary for eucariotic cells are supposed to af- 
fect mainly the motility and death of the cancer cells. So, 
the nutrients are divided into two groups: essential and 



and 



^ = V^TV - a^NcT,, ~ Xno'Nct, (4) 



^ = V^M - a'Ma^ ~ Xua^Ma, (5) 
ot 
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for the diffusion equations describing the nutrients con- 
centration fields. In addition, the boundary condition on 
the capillary vessel becomes N{x = 0) — M{x = 0) = 1 
and a value A = 1 is defined. 



ac{x) — 1 and CFd{x) = 1 when CTc vanishes. The cell death 
probability Pdei is determined by the concentration per 
cancer cells of the non-essential nutrients M present on 
the microenvironmcnt of the selected cell: 



C. Cell dynamics 

Each tumor cell can be selected at random, with equal 
probability, and carry out one of three actions: 

(1) division. Cancer cells divide by mitosis with prob- 
ability Pdiv If the chosen cell is inside the tumor, its 
daughter will pile up at that site, and adx) <Jc(x) + 1. 
Otherwise, if the selected cell is on the tumor border, its 
daughter cell will occupy at random one of their nearest 
neighbor sites x' containing a normal or a necrotic cell 
and, therefore, ac{x') = 1 and (Jn,d{x') — 0. The mitotic 
probability Pdiv is determined by the concentration per 
cancer cell of the essential nutrients N present on the 
microenvironmcnt of the selected cell: 



Pdei (x) = exp 



Pdiv{x) = 1 - exp 



N 



Cc tidiv 



(6) 



The Gaussian term is included in order to produce a sig- 
moid curve saturated to the unity, and the model param- 
eter 9div controls the shape of this sigmoid. 

(2) migration. Cancer cells migrate with probability 
Pmov A selected cell inside the tumor, at a site Xi, will 
move to a nearest neighbor site x' chosen at random. 
Thus, ac{x') ac{x') + \ and, clearly, (Tc(x) ac{x) — l. 
Otherwise, if the selected cell is on the tumor border, the 
invasion of a normal or necrotic nearest neighbor site will 
be dependent on the number of cancer cells present in 
the selected site. If in this site there is a single cancer 
cell, it migrates by interchanging its position with that 
of the invaded one. If there are other cancer cells in the 
same site of the one selected to move, the migrating cell 
will occupy the position of the normal or necrotic near- 
est neighbor cell, which, in turn, disappears. In terms of 
cell populations the migration of a cell on the tumor bor- 
der corresponds to the following operations: adx') = 1, 
ac{x) (Tc{x) - 1, Un^ix') = and an,d{.x) = 1 if 
ac{x) = 1. The probability of cell migration Pmov has the 
same functional form of Pdiv , but depends on the concen- 
tration of the non-essential nutrients M present on the 
microenvironmcnt of the selected cell and increases with 
the local population of cancer cells. So, 



Pmov [x) exp 



M 



(7) 



with the model parameter 9mov controlling the shape of 
this sigmoid. 

(3) cell death. Cancer cells die transforming in a 
necrotic cell with probability Pdei- Thus, (Jc{x) — > 



M 



(8) 



a Gaussian distribution whose variance depends on the 
model parameter 9dei- 

The cell dynamics rules used in our model take into ac- 
count that, as the cancer growth progress, cell migration 
increases near the border of the tumor due to the high 
availability of nutrients and the increase in the number 
of cancer cells, which release a series of enzymes (col- 
lagenases, metalloproteinases, etc.) responsible by the 
progressive destruction of the extracellular matrix. Also, 
in the regions where there is a high population density 
and an ineffective supply of nutrients via diffusion pro- 
cesses, the cell division is inhibited and, at the same time, 
the probability of cell death increases. But, under these 
rules cell growth and migration is possible even inside 
the tumor. Finally, the model parameters 9div^ dmov and 
9dei, which characterize the cancer cells response to nu- 
trient concentrations and embody complex genetic and 
metabolic processes, should be interpreted in terms of 
the underlying biochemistry and molecular biology, an 
still open problem. The other three model parameters 
a, Xn and Am, associated to the consumption of essen- 
tial and non-essential nutrients for cell proliferation by 
the normal and cancer cells, should be more easily deter- 
mined from biological experiments. 

It is worthwhile to notice that from the point of view 
of the so-called kinetic cellular theory, which provides a 
general framework for the statistical description of the 
population dynamics of interacting cells iQ, the local 
probabilities Pdivi Pmov and Pdei can be thought as an 
effective kinetic cellular model. 



D. Computer implementation 

The growth model simulations were implemented us- 
ing the following procedure. At each time step, the dif- 
fusion Eqs. (4) and (5) are numerically solved in the 
stationary state through relaxation methods, providing 
the nutrient concentration at any lattice site. Then, 
Nc'{t) cancer cells are sequantially selected at random 
with equal probability. For each one of them, a tentative 
action (division, death or movement) is chosen at ran- 
dom with equal probability and the time is incremented 
by At — \/Nc{t). The selected cell action will be im- 
plemented or not according to the correspondent local 
probabilities determined by Eqs. (6), (7) or (8). If the 
selected cell divides or die, therefore changing the num- 
ber of cancer cells which consume nutrients, we solve the 
diffusion equations in a small neighborhood of linear size 
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I = 20 centered in the altered site. This is done in or- 
der to take into account these local perturbations and 
to speed up the computer algorithm, since the number 
of numerical iterations need to solve the diffusion equa- 
tion is proportional to L^. At the end of this sequence of 
Nc{t) tentatives, a new time step begins and the entire 
procedure (solution of the diffusion equations and appli- 
cation of the cell dynamics) is iterated. The simulations 
stop if any tumor cell reaches the capillary vessel. In 
all the simulations, the exact solutions for the stationary 
diffusion equatios in the absence of tumor cells were used 
for the nutrient concentration fields at i = 0. 



III. SIMULATIONAL RESULTS 

In Figure 1 the most commonly observed morpholo- 
gies in tumor growth such as papillary, compact and 
disconnected are shown. Disconnected patterns, typical 
of round cells neoplasies such as lymphoma, mastocy- 
toma, and plasmacytoma, correspond to transient sim- 
ulated patterns in which cancer cells have high motility 
but a low mitotic rate. In turn, if the cell migration 
is very small and the high mitotic rate demands large 
amounts of essential nutrients, then the simulated pat- 
terns exhibit finger-like shapes similar to the papillary 
morphologies found in epithelial tumors, such as basal 
and skin cells carcinomas, and hepatomas. If nutrient 
availability in the tissue is further reduced by an increas- 
ing cell consumption, the papillary-like patterns become 
progressively thin and similar to the chords or filaments 
of cells which constitute one of the hallmarks of the tri- 
choblastoma morphology. Finally, under high nutrient 
supply, which means low cell consumption of both nu- 
trient types, the simulated patterns are compact such as 
those observed for solid tumors. 

Typical patterns generated by the present model are 
shown in Figure 2. As one can see, nonspherical mor- 
phologies growing towards the capillary vessel are ob- 
served, according to the rigorous results for moving 
boundary problems of cancer growth pO| . The nutri- 
ent consumption by normal and cancer cells, controlled 
by the model parameters a, Xn and Xm, plays a central 
role in morphology determination. For small values of 
these parameters, corresponding to growth conditions in 
which individual cells demand small nutrients supplies, 
the patterns tend to be compact and circular. However, 
if the mitotic rate of cancer cells is small due to the large 
amount of nutrients demanded for cell division, gener- 
ating a significant nutrient competition, these compact 
patterns progressively assume papillary-like morpholo- 
gies. At high nutrient consumption rates these papillary 
patterns become the rule and, for low cancer cell divi- 
sion, continuously transform in thin tips, filaments or 
chords of cells. Also, the smaller the a and Xn values 
are, the larger is the fraction of necrotic cancer cells for a 
fixed cellular response leading to cell death (controlled by 
Odei)- In particular, this fraction is smaller for papillary 



patterns than compact ones, suggesting that the optimal 
growth morphology under strong nutrients limitation is 
fractal. Finally, the higher cell migration is, more homo- 
geneous are the patterns, faster is the tumor growth, and 
smaller is the fraction of necrotic cells. 

The tumor patterns generated by the present model 
were characterized by its gyration radius Rg, total num- 
ber of cancer cells Nc-, and number of cells on tumor 
periphery S (including the surface of holes, if any). The 
gyration radius Rg is defined as 



(9) 



where n is the number of sites occupied by the pattern 
(cancer or necrotic cells) and is the distance of the 
occupied site i from the tumor mass center. These quan- 
tities could be related to clinically important criteria such 
as progress curves, rate of growth (volumetric doubling 
time) at given radii, proliferative and necrotic fractions of 
the tumor. In medicine, these data are used to determine 
tumor's malignancy and its prognosis. The obtained re- 
sults are summarized in Table 1. 

Despite of the detailed microscopic mechanisms of 
growth, for all the simulated patterns the progress in 
time of cancer cell populations follows a Gompertz curve 




Ncit) = Aexp[-exp(-^fc(i-t,))] 



(10) 



as one can see in Figure 3. Moreover, the tumor gyra- 
tion radius Rg as well as the number of peripheral cells 
S also exhibit a Gompertz growth, contrary to the lin- 
ear regimes observed for Rg and even S in our previous 
non nutrient-limited model Jl6t . So, the present model 
reveals that the Gompertz law of growth for the cancer 
cell populations and tumor size is a robust emergent fea- 
ture of cancer dynamics under nutrient competition. It is 
important to notice that linear or power law fittings for 
the growth in time of Rg and S can also be very satisfac- 
tory for several simulated patterns, therefore explaining 
our previous findings fl^ and the recent experimental 
observations of brain tumors grown in vitro |21 . How- 
ever, only the Gompertz law can provide good fittings 
for all the three quantities Nc, Rg and S. 

On the other hand, as a function of the total num- 
ber of cancer cells both Rg and S obey power law scal- 
ing given by Rg ^ and S ^ Nq, respectively. For 
solid patterns these exponents are v ~ 0.5 and a ~ 0.5, 
corresponding to effective circular and non fractal pat- 
terns. As the nutrient consumption increases the pat- 
terns tend to papillary-like shapes for which the exponent 
a increases towards the value 1 and the exponent v varies 
in the range [0.50,0.60], indicating a fractal morphology 
for the tumor. It is important to mention that increasing 
cell motility contributes to round and homogenize these 
patterns, progressively destroying their fractal features. 
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Also, the tumor growth is faster for higher cancer cell 
migration. 

An interesting result, shown in Figure 4, is the exis- 
tence of a necrotic core in the center of the simulated 
tumors for high nutrient consumption or cell division 
rates. As observed in real tumors and in vitro multi- 
cell spheroids , a simulated pattern consists of three 
distinct regions: a central necrotic core, an inner rim of 
quiescent cancer cells and a narrow outer shell of prolifer- 
ating cells. These different regions are evident in Figure 
5 where the cancer cell density and average cell division 
rate are plotted for a longitudinal cut across the growth 
pattern. As one can see, both the cancer cell density and 
mitotic rates have neat sharp maxima at the tumor bor- 
ders in front and opposed to the capillary vessel. Notice 
that the peaks for division rates are significantly more 
narrow than those for cancer cells density, demonstrat- 
ing that the proliferative fraction of the tumor comprises 
just a small part of the cancer cells localized at the tumor 
border. 

However, in the present nutrient-limited model the dis- 
connected patterns common in cancers of round cells 
(Fig. le) correspond to transient behaviors of the model 
for low mitotic rates and high cell migration of cancer 
cells. Also, the ramified morphology of a trichoblastoma 
seem in Fig. 6d could not be qualitatively reproduced by 
the model. A worthwhile feature of this pattern is the 
presence of "leaves" , the growth of each one clearly influ- 
encing the others. The simple introduction of cell motil- 
ity sensitive to the nutrients gradients failed to generate 
stationary disconnected or ramified patterns. Therefore, 
in addition to the nutrient field, it appears that a chemo- 
tactic interaction among cancer cells guiding their mi- 
gration must be in action. The nature of this biological 
interaction and the results of its simulation will be fo- 
cused in the next section. 



IV. GROWTH FACTORS 



^ = DnV^N - jNan ~ X^Nac , (11) 
ot 

in which 7 and A7 are the nutrients consumption rates of 
normal and cancer cells, respectively. The boundary con- 
ditions are the same described in Section 2. In turn, the 
growth factors (GFs) concentration obeys the diffusion 
equation: 



dG 
'dt 



^ DgV'^G ~k^G + VacN{GM -G) , (12) 



which includes the natural degradation of GFs, also im- 
posing a characteristic length ~ 1/fc for GFs diffusion, 
and a production term increasing linearly with the lo- 
cal nutrient concentration up to a saturation value Gj\/. 
Therefore, we are assuming that the release of GFs in- 
volves complex metabolic processes supported by nutri- 
ent consumption. The boundary conditions satisfied by 
the GFs concentration field is G{x, t) = at a large dis- 
tance (d > 2/fc) of the tumour border . 

Again, the number of parameters in equations (9) and 
(10) can be reduced by using the new dimensionless vari- 
ables 



Avi N_ G_ 



Using these new variables in Eqs. (9) and (10) and omit- 
ting the primes we obtain 



dN 
'dt 



(14) 



and 



From the biological point of view the chemotactic re- 
sponse to growth factors released by cancer cells seems to 
be, in addition to nutrient supply, another central feature 
in cancer development. The reciprocal influence among 
cancer cells mediated by autocrine and paracrine growth 
factors, motility factors, etc. influence the microenviron- 
ment of each cancer cell and hence its division and mi- 
gration. Thus, both diffusive fields (nutrients and growth 
factors) determine the local probabilities for cancer cells 
divide, migrate and die. 



A. The model 

In order to investigate the role of growth factors in a 
nutrient-limited growth, we simplify our model by con- 
sidering a single nutrient field described by the diffusion 
equation: 



dG 
'dt 



DV^G - k^G + TacN{l - G) 



(15) 



for the diffusion equations. The boundary condition for 
the nutrients on the capillary vessel becomes N{x = 0) = 
1 and a value A = 1 is defined. In addition, at the 
stationary state the parameter D in Eq. (13) can be 
put equal to the unity by rescaling the parameters k and 
r. So, the diffusion equations for the nutrients and GFs 
involve four parameters. 

On the other hand, the cell dynamics has essentially 
the same rules used in Section 2, but with different cell 
action probabilities. The single change introduced is that 
after cell division the daugther cell stands at the same 
site occupied by its mother. Since nutrients are essential 
to the large protein and DNA synthesis necessary to cell 
mitosis and GFs act are mitotic inductors, the proposed 
form to Prf,,, is: 
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Pdiv{x) = 1 - exp 



CTc 



(16) 



The parameter N* determines the nutrient-poor level be- 
low which the cancer cells reproduction is inhibited. 

Cell migration involves large citoskeleton reorganiza- 
tions which consume energy and is facilitated by GFs 
which destroy the extracellular matrix and the adhesivity 
structures between normal cells. Thus, Pmov is written 
as: 



-^mov(*^) *^ ) — 1 exp 



N{x)G{x)[G{x) -G{x')] 



On 



(17) 



implying that a cell migrates in a gradient-sensitive way 
towards sites where the GFs concentration is lower than 
that in its starting point. 

Finally, cell death is produced by the lack of nutrients: 



Pdei (x) = exp 



N 



CTc Odel 



(18) 



B. Results 

In Figure 6 typical compact, ramified and disconnected 
simulated patterns are shown. The ramified structure 
shown in Fig. 6b should be compared with the pattern 
of a trichoblastoma exhibited in Fig. 6d. In contrast to 
the compact fingers of the papillary patterns of the pre- 
vious section, in these ramified morphologies the tumor 
has fjords and tips similar to those abscrvcd in DLA pat- 
terns. We emphasize that without chemotactic signaling 
among cancer cells the nutrient-limited model can not 
generate stationary disconnected patterns. The reason is 
that in average GFs driven cell migration outwards, pro- 
moting the tumor expansion and, in consequence, gener- 
ating disconnected patterns for high cell motility. Again, 
the patterns were characterizcid by its gyration radius Rg, 
total number of cancer cells Nc ■ and number of cells on 
tumor periphery S. Essentially the same results were ob- 
tained for cancer progress in time, scaling relations and 
spatial structures exhibiting a central necrotic core, an 
inner rim of quiescent cells and a narrow outer shell of 
proliferating cells. 



V. CONCLUSIONS 

A nutrient-limited model for the growth of avascular 
tumors was investigated by numerical simulations. In 
its original version, cell proliferation, motility and death 



were locally regiilatcd by the concentration of nutrients 
supplied by a distant capillary vessel. These nutrients 
were divided into two groups, the first one associated to 
the usual metabolic cell needs, and the second essential 
to the synthesis of proteins and nucleic acids involved in 
cell division. The nutrients concentration fields were de- 
termined by solving the diffusion equation on the square 
lattice modeling the primary tissue. Our simulation re- 
sults show that the progress in time of the total number of 
cancer cells, tumor gyration radius and number of cells 
on the tumor border is described by Gompertz curves. 
The generated compact and papillary or fingcr-likc mor- 
phologies obey different scaling laws for the number of 

1 /2 

peripheral cancer cells. For compact patterns S ~ Nq 
as in the Eden model, whereas for papillary patterns the 
exponent in the power-law increases towards unity as the 
nutrient consumption increases, indicating a fractal mor- 
phology for the tumor. Since in this model version the 
cell migration is not driven by chemotactic signals se- 
creted by the cancer cells, cell motility contributes to 
round and homogenize the growth patterns. Also, the 
simulated tumors incorporate a spatial structure com- 
posed of a central necrotic core, an inner rim of cjuics- 
cent cells and a narrow outer shell of proliferating cells 
in agreement with biological data. 

In order to simulate disconnected and ramified tumor 
patterns, typical of round cells tumors and trichoblas- 
toma, a chemotactic interaction among cancer cells me- 
diated by growth factors was added to the competition 
for nutrients. Again, similar results were obtained for 
cancer progress in time and scaling relations. Thus, the 
Gompertz law emerges as a robust feature of the nutrient- 
limited model of cancer growth. 

Despite the encouraging results for the progress curves, 
scaling laws and growth patterns morphologies, includ- 
ing nonspherical symmetries, the major feature of the 
present work is an attempt in connecting the macroscopic 
diffusion equations for nutrients and/or growth factors 
to cell response and interactions at the microscopic scal- 
ing through an effective kinetic cellular model. Indeed, 
the local probabilities Pdiv, Pmov and Pdei describe in a 
stochastic way the dynamical processes occuring in cell 
populations as a response to the nutrient and growth fac- 
tors diffusive fields. Finally, further studies on angiogen- 
csis, therapy and tumor-host interactions using variants 
of the present model are under progress. 
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TABLE I: Morphology, progress curves and scaling laws for 

the patterns generated by the nutrient-limited model 
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FIG. 1: Commom morphologies observed in cancer growth, (a) Papillaxy pattern of a scamous papyloma, (c) a compact sohd 
basocellular carcinoma, (e) a disconnected pattern of a plasmacytoma, and (g) characteristic cell filaments of a trichoblastoma. 
All these histological patterns were obtained from dogs. The corresponding simulated patterns are shown in (b), (d), (f) and 
(h), respectively. 
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FIG. 2: Simulational results of the nutrient-limited cancer growth model. The patterns are organized as function of nutrient 
consumption rate a for normal cells and the multiplicative factor Ajv to the consumption rate of mitotic essential nutrients 
by the cancer cells. The remaining four parameters of the model were fixed in Am = 10, duv = 0.3, 9mov = oo (absence of 
cell migration) and Odei = 0.01. The patterns are drawn in a gray scale where the darker regions represent higher cancer cell 
populations. The tissue size is 500 x 500, with the initial "cancer seed" distant 300 sites from the capillary. The total number of 
cancer cells depends on tumor morphology and attains up to 2 x 10 ' for compact patterns. The simulated patterns arc compact 
for low Ajv values and become papillary or finger-like for high Ajv. For the same Ajv the patterns are more papillary for higher 
a. Since the capillary vessel provides a fixed nutrients supply, the dimensionless consumption rate of the normal tissue a set up 
the levels of available resources for which cancer cells compete. So, high a and/or Xn values correspond to the limit of strong 
nutrient competition. 
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FIG. 3: Gomportz growth of cancer cell population. The model parameters were fixed in q = 2/L, Am = 10, Xit 
Odiv = 0.3, 9mov = 2 and 6dei = 0.03. The solid line correspond to a nonlinear fitting with = 0.9995. 



200, 




FIG. 4: Simulated growth patterns exhibiting a necrotic core (in black). The model parameters are a = 2/L, Am = 25, 
6div = 0.3, 9mov = oo (without cell migration) and Odei = 0.03. The fraction of necrotic cells is smaller for papillary patterns 
than compact ones, suggesting that the optimal growth morphology under strong nutrients limitation is fractal. 
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FIG. 5: (a) Density of cancer cells Uc and (b) division rates Wdiv as a function of the distance from the capillary vessel along 
a tumor longitudinal cut. The model parameters were fixed in a = 2/L, Am = 25, Ajv — 50, 6div = 0.3, 6mov = oo (without cell 
migration) and ddei = 0.03. Sharp mgixima at the tumor borders near and opposed to the capillary vessel at x = are evident. 
Therefore, the proliferative fraction of cancer cells is distributed in a thin shell on the tumor border. 




FIG. 6: Simulated patterns of the nutrient-limited cancer growth model including the influence among cancer cells mediated by 
growth factors, (a) compact {k = 0.1 and 6mov ~ 1), (b) discomiected (fc — 0.025 and 6,nov = 0.1), and (c) ramified (fc = 0.025 
and 9mov = 1) morphologies. The remaining parameters of the model were fixed in a = 3/L, 9div = 0.5, A'"* = 9dei = 0.01, 
A = 5 and F = 10. The tissue size is 500 x 500, with the initial "cancer seed" distant 300 sites from the capillary and the total 
number of cancer cells is 5 x 10*. For comparison, a real ramified pattern observed in trichoblastoma is shown in (d). 



